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Starting with an expression, due originally to Einstein, for the shear vis- 
cosity r]{5(f)) of a liquid having a small fraction 5(f) by volume of solid partic- 
ulate matter suspended in it at random, we derive an effective-medium vis- 
cosity ri{(j)) for arbitrary cj) which is precisely of the Vogel-Fulcher form. An 
essential point of the derivation is the incorporation of the excluded-volume 
effect at each turn of the iteration 0„+i=0„-/-50. The model is frankly me- 
chanical, but applicable directly to soft matter like a dense suspension of 
microspheres in a liquid as function of the number density. Extension to a 
glass forming supercooled liquid is plausible inasmuch as the latter may be 
modelled statistically as a mixture of rigid, solid-like regions (0) and floppy, 
liquid-like regions (1-0), for increasing monotonically with supercooling. 

PACS number(s): 64.70.Pf, 45.70.-n, 66.20.+d, 83.80.Hj 

Extreme slow dynamics defines approach to the glassy state. At the 
macroscopic scale, it manifests as a rise of shear viscosity, typically by 15 
orders of magnitude, as that state is reached through supercooling of the 
glass forming liquid. The Vogel-Fulcher (VF) law describes that growth of 
viscosity [Ij.This work derives the VF law [2]. 

A striking feature of the VF law is the essential singularity, rather than a 
power-law divergence, of the shear viscosity at a temperature Tq. The relax- 
ation times, however, exceed the experimental time scale at what is identified 
as the glass transition temperature Tg > Tq, making thus the glass transition 
a kinetic crossover. This inverse exponential VF law is well known to hold 
for the fragile structural-glass forming liquids [1]. But, very significantly it 
is also obeyed by a broad class of soft-matter systems that exhibit the ex- 
treme slow dynamics [1]. This includes the purely mechanical systems, e.g., of 
weakly perturbed granular aggregates, where the degree of compaction and 
the perturbation strength, rather than mass density and temperature, are 
the relevant variables and the control parameter, and the underlying physics 
is that of jamming, or blocking, by rigid granular contacts [3-6]. And, simi- 
larly for the case of a dense suspension of microspheres [1]. Motivated by its 
ubiquity and universality, we have attempted a derivation of the VF law for 
a fluid-mechanical model of a liquid containing a volume fraction of solid 
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particulate matter suspended in it at random. It is an effective medium tlie- 
ory (EMT) along the line of Bruggemann's asymmetric EMT [7], but it goes 
beyond the mean field by incorporating the solid-solid exclusion explicitly in 
real space, which indeed is the essential point of our derivation. This, frankly 
mechanical model can, however, be re-interpreted as a model for the glass 
forming supercooled liquid inasmuch as the latter may be re-approximated 
as a statistical mixture of short-ranged rigidity (solid-like fraction (p) and the 
floppy liquid-like fraction (1-0). In our view the present work complements 
the other derivations of the VF law which are based on the idea of marginal 
scaling [6] and some simple exclusion models [8-11]. 

We start with the expression, due originally to Einstein [12], for the shear 
viscosity f]{S(f)) of a liquid containing a small volume fraction 50 of solid 
particulate matter suspended in it at random: 



where a, of order unity, is a fluid-dynamic dimensionless parameter specifying 
the particle shape and the flow boundary condition, and 50 = {4'^/3)a^6n 
assuming spherical particles of radius 'a' with 6n the number density. The 
physical basis of Eq. ((T)) is that in the steady state the rigid parts of the liquid 
move practically as complete wholes, and hence the effect of their existence 
is to diminish the thickness of the layer, through which momentum has to be 
transported by the mobile molecules, and thus to increase the viscosity [13]. 
We can iterate Eq. (Q) to a higher volume fraction 0, in the spirit of an 
EMT, by the recursion relation 



where the factor 1/(1-0) in the denominator on the right-hand side ensures 
that the elemental increment 6(f) is reckoned relative to the liquid-like volume 
fraction (1-0) remaining at the current stage of iteration. Now, proceeding 
to the limit 50 ^0, we obtain the differential equation 



r/(0) (1 + a50) 



(1) 



r/(0 + 50) = 7^(0) (1 + a 



(2) 




(3) 



with the solution 



r/(0)=r7(O)(l-0) 



(4) 
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that gives a power-law divergence for the effective shear viscosity ri{(j)). Here 
77(0) is the 'bare' viscosity of the pure hquid with = 0. Such a power-law 
temperature dependence is well known to follow from the viscosity feedback 
mechanism giving the Batchinski-Hildebrand law [14] (with a=l), or from 
the Mode Coupling Theory (MCT) [15] giving the critical behaviour (with 
a=2). Both these exponent values lie in the range for the parameter a as 
described below. 

Equation (0} giving this critical behaviour is, however, in error in that 
it mathematically ignores the physically important excluded- volume effect. 
The point is that the liquid fraction (1-0) in the denominator in Eq. Q must 
be replaced by the liquid fraction (1-0) as weighted by the probability that 
the incremental solid fraction 6(f), added at random, lands in it [16]. This, 
therefore, effectively replaces (1-0) by (1-0)^' Equation (jlj) then gets modified 
accordingly to 



drj a 

giving 



(5) 



77(0) = (r/(0) e"") e"/^^-'^). (6) 

The expression in Eq. © is already of the VP form as an inverse ex- 
ponential function of diverging essentially at =1. This, however, needs 
a refinement as dictated by the physics of the problem, namely, that the 
solid volume fraction need approach only the rigidity percolation thresh- 
old 00 (jl) in order to attain the three-dimensional rigidity. Therefore, (1-0) 
above must be displaced to (0o - 0). Thus, we finally have 

r^(0) = (r/oe-°/'^«) e"/(^«-'^\ (7) 

which tends to ri(0) for — >0 (pure liquid), and diverges as ^ 0o from 
below (the glassy state). 

Equation (|7j) is our main result. For the simplest case of spherical, non- 
spinning particles, we have [12,13] a = 2.5, while for particles free to spin, 
a=l. Also, we can estimate the rigidity percolation threshold [17] (0o'^)in 
three dimensions from its 2D value 0o^ ~ 0.80 by use of the simple relation 
03^ = 4/37ri/2(02D)3/2. We get 0^^ ~ 0.54. In Fig. (HJ, we have plotted 
77(0) against for the values of the parameters, a = 2.5 and 00^^= 0.54. This 
is essentially a universal curve. 
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While Eq. ((7j) is expected to be directly applicable to, e.g., a suspension 
of microspheres in a viscous liquid, its extension to the glass forming su- 
percooled liquids is plausible as indicated earlier. Then cj) must be regarded 
as a function of temperature, increasing monotonically as the temperature 
decreases. This will turn Eq. © explicitly into the VF form, or its variant, 
the Vogel-Tammann-Fulcher (VTF) law, rj{T) = rjo exp{DTQ/(T — Tq)) as 
T-^ Tq from above. 

We would now like to conclude with the following remarks. The above 
fluid-mechanical model physically implies that our derivation may apply more 
readily to the fragile rather than to the strong (network forming) liquids. As 
noted above, the numerical value of a occurring in Eq. ((7j) depends on the 
particle shape [12] (taken to be spherical here), and on whether the particles 
are free to spin [a = 1) or not [a = 2.5) in the presence of a shear rate. 
This can make the parameter a temperature dependent, with the higher 
value a = 2.5 appropriate to the lower temperatures. With the solid-like 
volume fraction now become a function of temperature, and, therefore, a 
thermodynamic parameter, the equation ((7j) shows how the shear viscosity (a 
transport property) is actually controlled by the thermodynamics: The ther- 
modynamically controlled liquid-like fraction (0o -0) acts as an idler taking 
up the shear rate. This is the simplest realization of a viscosity amplification 
that underlies the macroscopic slow dynamics described by the Vogel-Fulcher 
law, where the idling liquid-like fraction essentially retains its 'bare' low value 
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FIG. 1. Plot of normalized shear viscosity T]{(f))/r](0) against the solid-like 
volume fraction from Eq. (|7j) derived in the text, for a = 2.5. Here 0o is 
the rigidity percolation threshold, and (j)g marks the point {(j)g)/ri(0)= 10^^. 
The regime 0/0 /0g is nominally the supercooled liquid; 0^/0 i4>(0) the glassy 
liquid; and 0^0o the rigid glassy solid. 
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